Exponential speed-up with a single bit of quantum information: 
Testing the quantum butterfly effect 
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We present an efficient quantum algorithm to measure the average fidelity decay of a quantum 
map under perturbation using a single bit of quantum information. Our algorithm scales only 
as the complexity of the map under investigation, so for those maps admitting an efficient gate 
decomposition, it provides an exponential speed up over known classical procedures. Fidelity decay 
is important in the study of complex dynamical systems, where it is conjectured to be a signature 
of quantum chaos. Our result also illustrates the role of chaos in the process of decoherence. 
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"The flap of a butterfly's wings in Brazil can set off 
a tornado in Texas" . This butterfly effect illustrates the 
canonical feature of chaotic systems: they display ex- 
treme sensitivity to their initial conditions. In a chaotic 
regime, phase space trajectories diverge exponentially in 
time at a rate governed by the largest Lyapunov exponent 
of the system. Isolated quantum systems cannot display 
the butterfly effect, since unitary evolution preserves dis- 
tances between states. Extensive research over the past 
two decades has been devoted to examining other man- 
ifestations which can be used to distinguish the regular 
and chaotic regimes of quantum systems. While many 
signatures of quantum chaos have been proposed, their 
validity relies mostly on vast accumulations of numeri- 
cal evidences. Furthermore, obtaining conclusive results 
using these measures requires manipulating data whose 
size scales as the dimension (N) of the system's Hilbcrt 
space - that is, exponentially with the number of qubits 
K required to simulate the system. In this article, we 
demonstrate a quantum algorithm to evaluate one such 
signature — the average fidelity decay — with a single bit 
of quantum information, in a time that scales as poly(A'). 

Fidelity decay was initially proposed as a signature of 
chaos by Peres |lj , and has since been extensively investi- 
gated 0, 0] ■ The closest quantum analogue to the (purely 
classical) butterfly effect, fidelity decay measures the rate 
at which identical initial states diverge when subjected to 
slightly different dynamics. The discrete time evolution 
of a closed quantum system can be specified by a unitary 
operator U, where p(r n ) = U n po(W) n . To examine fi- 
delity decay, we construct a slightly perturbed map U p , 
where U p = UP with P = cxp{— iSV} for some small 
S and a hermitian matrix V. It is conjectured that the 
overlap (or fidelity) 
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(as a function of the discrete time n) for regular and 
chaotic dynamics: chaotic dynamics will display expo- 
nential fidelity decay, while regular dynamics will pro- 
duce polynomial fidelity decay. Actual results to date 
show behavior rather more complex than the preceding 
simple conjecture. There are various regimes governing 
the decay rate; fidelity decay depends on the perturba- 
tion strength S, and also on the degree of correlation 
between the eigenbasis of U and the eigenbasis of the 
perturbation P Fidelity decay remains a powerful 
diagnostic of chaotic behavior, but calculating it is com- 
putationally hard. Furthermore, because F n (ip) gener- 
ally shows large fluctuations over time, it is in practice 
necessary to average F n (tp) over a random set of initial 
states if) to determine its decay rate, thus increasing the 
numerical burden. 

Several classically hard problems can be solved in poly- 
nomial (in A) time on a quantum computer. Since fully 
controllable and scalable quantum computers are still 
quite a ways in the future, algorithms which can be per- 
formed on a less-ambitious quantum information proces- 
sor (QIP) are of great interest. A QIP is a quantum 
device which may fail to satisfy one or more of DiVin- 
cenzo's five criteria, but can nonetheless carry out in- 
teresting computations 0. Of particular interest to us 
is deterministic quantum computation with a single bit 
(DQC1) 0, a model of quantum information processing 
which is believed to be less powerful than universal quan- 
tum computation and which is naturally implemented by 
a high temperature NMR QIP Q. In this model, uni- 
versal control over all qubits is still assumed, but state 
preparation and read-out are limited. The initial state of 
the (A + l)-qubit register is 
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between initially identical states ip undergoing slightly 
different evolutions, U and U p , should decay differently 



i.e., the first qubit (called the probe qubit for rea- 
sons which will become clear) is in a pseudo-pure state, 
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whereas the other K qubits are in the maximally mixed 
state. Furthermore, the result of the computation is ob- 
tained as the noisy expectation value of a z on the probe 
qubit. The variance of a z is determined by i) the po- 
larization 7 of Eq. [21 (independent of the size of the reg- 
ister) and ii) the inherent noise of the measuring pro- 
cess. Hence, (a z ) can be estimated to within arbitrary 
e with a probability of error at mostn by repeating the 
computation 0(log(l/p)/e 2 ) times The value of 7 
in high-temperature NMR is independent of the size of 
the register because only a single qubit needs to be in a 
pseudo-pure state. The "inherent noise" receives contri- 
bution from both electronic noise and statistical fluctua- 
tions due to the finite sample size. 

While it has been known for some time that the dy- 
namics of some quantized chaotic systems can be ef- 
ficiently (i.e^, in poly(A) time) simulated on quantum 
computers ||| , it was shown only recently [J, [lfj, |ll| that 
this ability can also be used to efficiently evaluate cer- 
tain proposed signatures of quantum chaos. In Ref. 1 1 Ol ] , 
an efficient quantum circuit is constructed to evaluate 
the coarse grained local density of states (LDOS) — the 
average profile of the eigenstates of U over the eigenba- 
sis of U p — which is believed to be a valid indicator of 
chaos and is formally related to fidelity decay via Fourier 
transform j^. In Ref. fl, an efficient procedure to es- 
timate the fidelity decay using the standard model of 
quantum computation is presented. Finally, in Ref. [ill ], 
a DQC1 circuit is presented to estimate the form factors 
t n = \Tr{U n }\ 2 of a unitary map U which, under the 
random matrix conjecture (see [12j | and refs. therein), is 
a good signature of quantum chaos. The proposed algo- 
rithm offers only a quadratic speedup, but since entangle- 
ment is very limited in DQC1 |l3|, this result raises doubt 
about the common belief that massive entanglement is 
responsible for quantum computational speed-up [Til ]. 

Drawing upon all this previous work, we will now con- 
struct an efficient DQC1 algorithm to evaluate the av- 
erage fidelity decay associated with any pair of unitary 
operators U and U p , provided thev can be implemented 
efficiently, e.g. as those of Refs. p. We begin by prov- 
ing a crucial identity required to implement the efficient 
algorithm. 

Let f(ip) be a complex-valued function on the space 
of pure state of a A-dimcnsional quantum system. We 
denote its average by f(ip) = J f(ip)dip, where dip is the 
uniform measure induced by the Haar measure, such that 
J dip = 1. For sake of compactness let (ip\A\ip) = (A)^. 
Theorem: Let A, B, C, ... be £ linear operators on a TV- 
dimensional Hilbert space. Then 



Tr 



(£) 

where P s is the projector on the symmetric subspace of 
£ systems, see Ref. for details on P^p . 
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Proof: First, note that 

(A)^(B)^(C)^ . . . = Tr { \iP)(iP\® e (A ® B ® C . . .)} . 
Therefore, the average over the pure states ip yields, 
{|V)W>|® £ (A(g)B®C*...)}. 



Tr 



Since \ip){ip\® e annihilates any state which is antisym- 
metric under interchange of two of the £ systems, and 
is by construction symmetric under such interchange, it 
must be proportional to the projector Pg onto the sym- 
metric subspace. To establish the theorem it is suffi- 
cient to find the proportionality factor A between these 
two quantities. Letting A = B = C= ... = 'R i we 
get 1 = Tr^W^} = ATr{pf } = A( JV+ /" 1 ) (see 
Ref. [lij)j which completes the proof. □ 
A useful corollary to this Theorem for any specific I 
can be obtained by expanding Pp^ in Eq. [21 In the case 
£ = 2, it reads 
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Similar expressions can be derived for £ > 2, which in- 
volves the properly normalized sum of all combinations 
of traces of products and products of traces. 

To arrive at our algorithm, it is sufficient to write the 
average fidelity as F n (ip) = ({U n )m; l )^((U^U n }^, and 
apply the identity from Eq. 0]to obtain 



|Tr{([/ n )t£/™}| 2 +N 
N 2 + N ' 



(5) 



The specific form of our theorem with £ ~ 2, unitary 
A, and B = A^ was discovered by M., P., and R. 
Horodecki , but our proof simplifies the presentation. 
An efficient DQC1 algorithm to evaluate the trace of any 
unitary operator [here, (U n )' ! Up], provided that it ad- 
mits an efficient gate decomposition, was presented in 
Ref. 0| ■ If the perturbed map takes the form U p = UP 
for some unitary operator P (e.g., P ~ exp{— iSV} as 
above), the circuit can be further simplified into the one 
illustrated on Fig. ^ 

We now analyze the complexity of our algorithm. We 
assume that U and U p admit e-accurate gate decom- 
positions whose sizes grow as L(K,e) £ poly(A^, 1/e). 
This implies that the controlled version of these gates 
also scale as L(K,e) 0|. We see from Eq. [5] that 
the variance of F n (ip) is at most twice the variance of 
Tr{{U n yU%}/N. Therefore, the overall algorithm - esti- 
mating F n (ip), to within e, with error probability at most 
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FIG. 1: Quantum circuit evaluating the average fidelity 
F n (ip) between the perturbed and unperturbed maps U and 
Up = UP. The gates Rk are ir/2 rotation in the Bloch sphere 
around axis k = x or y. When k is set to i, we get the real 
part of Tr{(U n )^Up}/N while k = y yields the imaginary 
part. The unitary operator P is applied conditionally: when 
the probe qubit is in state |1), the unitary P is applied to the 
lower register while no transformation is performed when the 
state of the probe qubit is |0). 



p - requires resources growing as L(K, e)n log(l/p) /e 2 , so 
it is efficient. (The range in n over which the decay is 
studied should be independent of the system's size.) This 
algorithm thus provides an exponential speed-up over all 
known classical procedures and uses a single bit of quan- 
tum information. Furthermore, it eliminates any cost of 
averaging the fidelity over a random set of initial states, 
as this averaging is done directly. 

In order to implement certain unitary maps on K 
qubits efficiently, it is necessary to introduce a number 
K a of ancillary qubits (a "quantum work-pad" ) in the 
fiducial state IV'o}- Ancillary qubits in pseudo-pure states 
can be used in the DQC1 setting. As a first step of the 
computation, part of the polarization of the probe qubit 
of Eq. can be transferred to ancillas initially in maxi- 
mally mixed states. Thus, as long as the size K a of the 
work-pad is at most poly-logarithmic in K . the algorithm 
remains efficient. 

Perhaps the most surprising feature of the quantum al- 
gorithm as it is presented in Fig.^is that the probe never 
gets entangled with the system throughout the compu- 
tation. To show this, consider a generalized version of 
the circuit of Fig. ^where the P's and the £/'s are free 
to differ at each iteration, i.e. at step j, we apply Pj 
conditionally on the probe qubit, followed by Uj. This 
generalization is necessary since the controlled P gate 
will in general be decomposed as a sequence of elemen- 
tary controlled and regular gates . Initially, the probe 
qubit is in state a\0) + f3\l). After k steps, the state of 
the QIP is 



! |0X0| ® 11 + a/3*|0Xl| ® S + 



+ a*/3|l)<0|®S+|/3| 2 |l)<l|®ll} 
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where S = U k P k . . . U 2 P 2 U 1 P 1 Ulul . . . U\. Decompos- 
ing this state in the eigenbasis of the unitary matrix 

S\(j>j) = e is i\<f>j), we get 



Pk 



-^5Z|ayX a jl®l&X0j 



(T) 



where |ay) = a|0) + /3e tSj \l); the state is separable. Its 
separability supports the point of view that the power of 
quantum computing derives not from the special features 
of quantum states — such as entanglement — but rather 
from fundamentally quantum operations |l9l |. 

Our algorithm also illustrates why chaotic environ- 
ments are expected to produce decoherence more rapidly 
than integrable ones |2(J. Consider the probe qubit of 
Fig. n as a quantum system interacting with a complex 
environment consisting of K two-level systems. After 
a "time" n, the state of the system is given by tracing 
out the K environmental qubits from Eq. The diago- 
nal elements of the reduced density matrix |cv| 2 and |/3| 2 
are left intact while the off-diagonal elements a(3* and 
a* (3 are decreased by a factor | Tr { 5" } | which is roughly 



equal to y F n {ip). Thus, for an environment with chaotic 
dynamics, the system will decohere at an exponential 
rate, whereas the rate of decoherence should be slower 
for non-chaotic environments. This analogy also provides 
a very simple example of decoherence without entangle- 
ment |2lj . 

On the circuit of Fig. ^ only the perturbation gates 
P are conditioned on the state of the probe qubit. This 
suggests a dual interpretation of the algorithm as quan- 
tum circuit and quantum probe. On the one hand, U 
could be a known unitary transformation which is be- 
ing simulated on the lower Ji'-qubit register over which 
we have universal control. Then, the gate U would sim- 
ply be decomposed as a sequence of elementary gates as 
prescribed in Refs. || for example. On the other hand, 
the lower register could be a real quantum system under- 
going its natural evolution U which might not even be 
known. Then, the probe qubit should really be regarded 
as a probe which is initialized in a quantum superposi- 
tion, used to conditionally kick the system, and finally 
measured to extract information about the system under 
study. In this case, it is not necessary to have universal 
control over the lower register (the quantum system), we 
must simply be able to apply a conditional small unitary 
transformation to it. 

Finally, Eq. [S] provides a useful numerical tool that 
can be used to compute the exact average fidelity in- 
stead of estimating it by averaging over a finite ran- 
dom sample of initial states. In Ref. Q, fidelity de- 
cay was illustrated on the quantum kicked top map 
Uqkt = exp{— iir J y /2} exp{— ik J^/j} acting on the 
N = 2j + 1 dimensional Hilbert space of angular mo- 
mentum operator J. The chosen perturbation operator 
was P = Ylf = i exp{— i5o- 3 z /2}, a collective rotation of all 
K qubits of the QIP by an angle 5. The decay rate (gov- 
erned by the Fermi golden rule in this regime) was eval- 
uated to T = 2.50(5 2 for this perturbation 0. Fni^*) was 
estimated in both chaotic (k = 12) and regular (k = 1) 
regimes of the kicked top by averaging over 50 initial 
states. We reproduce these results on Fig.|5|and compare 
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them with the exact average Eq. [5] and theoretical pre- 
diction e~ r ™. The random sample is in good agreement 
with the exact average except that the former shows fluc- 
tuations. Furthermore, the decay in the chaotic regime is 
in excellent agreement with the Fermi golden rule. While 
the decay is slower in the regular regime, it is not clear 
from these results that it is not exponential. 
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FIG. 2: Fidelity decay F n (ip) averaged over 50 initial com- 
putational basis states for Uqkt in a regular regime (k = 1, 
squares) and chaotic regime (k = 12, circles). The dashed 
lines represent the exact average Eq. |S]and the full line shows 
the exponential decay at the Fermi golden rule rate P. 

We have presented an efficient quantum algorithm 
which computes the average fidelity decay of a quantum 
map under perturbation using a single bit of quantum 
information. The quantum circuit for this algorithm es- 
tablishes a link between decoherence by a chaotic envi- 
ronment and fidelity decay. Using a special case of our 
theorem, we have numerically evaluated the exact aver- 
age fidelity decay for the quantum kicked top, and found 
good agreement with previous estimations using random 
samples. Although we have mainly motivated our algo- 
rithm for the study of quantum chaos, we believe that it 
has many other applications such as characterizing noisy 
quantum channels and computing correlation functions 
for many-body systems. We have also shown that our 
algorithm can be viewed as a special experiment where a 
quantum probe is initialized in a superposition and used 
to conditionally kick the system under study. This type 
of quantum information science byproduct might open 
the horizon to new types of experimental measurements 
where a small QIP is used to extract information from 
the quantum system under study. Finally, the effective 
speed-up despite the limited presence of entanglement — 
in particular its complete absence between the quantum 
probe and the mixed register — is a step forward in our 
understanding of the origin of quantum-computational 
speed-up. 

We thank J. Emerson, G. Milburn, J. P. Paz, and W.H. 
Zurek for helpful discussions. We also acknowledge the 
Benasque Center for Science where this work was ini- 



[3 
[4 
[5 
[6 

[' 

[« 
[9 

[10 

[11 
[12 

[13 
[14 

[15 

[10 

[17 
[18 

[19 

[20 
[21 



A. Peres, Phys. Rev. A 30 1610 (1984). 

R.A. Jalabert and H.M. Patawski, Phys. Rev. Lett 

86 2490 (2001); F. Cucchietti, C.H. Lewenkopf, E.R. 

Mucci olo, H. Patawski, and R.O. Vallejos, arXiv: 

|nlin.CD/0111051| (2001); G. Benenti and G. Casati, 

arXiv: jquaiiTp h/0112060*] (2001); T. Prosen and M. 

Znidaric, J. Phys. A 35 1455 (2002); 

P. Jacquod, PG. Silvestrov, and C.W.J. Beenakker, 

Phys. Rev. E 64 55203 (2001). 

J. Emerson, Y. Weinstein, S. Lloyd, and D. Cory, Phys. 
Rev. Lett. 89 284102 (2002). 

R. Blume-Kohout, CM. Caves, and I.H. Deutsch, Foun- 
dations of Physics 32, 1641 (2002) 

E. Knill and R. Laflamme, Phys. Rev. Lett 81 5672 
(1998). 

D.G Cory et al., Fortschr. Phys. 48, 875 (2000). 
P. J. Huber, Robust statistics (Whiley, New York, 1981). 
R. Schack, Phys. Rev. A 57, 1634 (1998); B. Georgeot 
and D.L. Shepelyansky, Phys. Rev. Lett. 86, 2890 (2001); 
G. Benenti, G. Casati, S. Montangero, and D.L. Shep- 
elyansky, Phys. Rev. Lett 87, 227901 (2001). 
J. Emerson, S. Lloyd, D. Poulin, and D. Cory, arXiv: 
|quant-ph/0308 164 (2003). 

D. Poulin, R. Laflamme, G.J. Milburn, and J. P. Paz, 
Phys. Rev. A 68 022302 (2003). 

F. Haake, Quantum Signatures of Chaos, Spriger-Verlag, 
Berlin (2001). 

K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 
stein, Phys. Rev. A 58 883 (1998); S.L. Braunstein, CM. 
Caves, R. Jozsa, N. Linden, S. Popescu, and R. Schack, 
Phys. Rev. Lett. 83 1054 (1999). 

A. Ekert and R. Jozsa, Phil. Trans. R. Soc. Lond. A 1769 

(1998); Jozsa, R. in The Geometric Universe edited by 

S. Huggett, L. Mason, K. P. Tod, S. T. Tsou and N. 

Woodhouse, Oxford University Press (1998). 

A. Barenco, A. Berthiaume, D. Deutsch, A. Ekert, R. 

Jozsa, and C. Macchiavello, SIAM J. Comput. 26 1541 

(1997). 

M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 
Rev. A 60 1888 (1999). 

C. Miquel, J. P. Paz, M. Saraceno, E. Knill, R. Laflamme, 

and C. Negrevergne, Nature 418, 59 (2002). 

A. Barenco, C.H. Bennett, R. Cleve, D.P. DiVincenzo, 

N. Margolus, P.W. Shor, T. Sleator, J.A. Smolin, and H. 

Weinfurter, Phys. Rev. A 52, 3457 (1995). 

R. Laflamme, D.G. Cory, C. Negrevergne, and L. Viola, 

Quant. Info, and Comp. 2 (2001); D. Poulin, Phys. Rev. 

A 65 42319 (2002). 

R. Blume-Kohout and W. H. Zurek, arXiv: 
|quant-ph/0212153| (2002) . 

J. Eisert and M.B. Plenio, Phys. Rev. Lett. 89 137902 
(2002). 



